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Abstract 


The objective of the present study was to obtain an 
empirical correlation to relate the ultimate depth to which 
a high pressure water jet can drill or pierce a hole into an 
Oil sands formation to the jet parameters. To this end a 
laboratory apparatus was constructed to study the 
characteristics of confined, continuous, high pressure water 
jets in a down-hole type geometry. Jet dynamic pressures up 
toesSeMPacdnd «nozzle srzessofs3<2pe2 04, cand 1.6 mm-in 
diameter were used. Horizontal tubes of 102, 76, 51 and 38 
mm diameter and a porosity of 5% for the 51 mm tube diameter 
were used in the tests to simulate the hole cut by the jet. 
Standoff distances ranged from 0.3 to 1.7 meters. 

Flow visualization observations and records of the 
pressure signal at the bottom of tube revealed that the 
water jet produced discrete impacts. From this observation 
it was concluded that it 1S most appropriate to model the 
jet drilling as a stochastic process and to analyze the 
pressure Signal statistically. A correlation to predict the 
ultimate depth of penetration was developed based on the 
probability of the pressure at the cutting face exceeding a 
Critical value for non-negligible cutting rates to occur. A 
linear relationship between the cutting rate and the above 
probability was assumed and a second correlation predicting 


the depth of penetration with time was developed. 
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The predicted depth Beipeneereicn agreed closely with 
the field test data for jetting times up to 60 seconds. No 
data 1S available for jetting times beyond 60 seconds. The 
model predicts a lower value than the correlation derived 
from experimental work usSing recompacted oil sands (Chau, 
1980). 

The model could not be effectively evaluated due to the 
limited field test results. Modelling of the jet cutting as 
a stochastic process appears to be an appropriate approach 
and 1S considered to give more valuable information than any 


cutting models based on the average or peak pressure. 
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nozzle diameter 


tube diameter 


jet dynamic pressure at the nozzle 


jetting time 

time to reach ultimate depth 

average velocity of jet at the nozzle 
depth of penetration 

ultimate depth of penetration 

density of water 

volume flow rate 


pressure at the cutting face 


Critical pressure at 011 sands surface below 


which no further penetration occurs 


Inetital “curting) Tate 


probability that the pressure p exceeds Pc 


experimental constants 


cutting rate 
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1. Introduction 


1.1 Background Information. 

Oil sands, also known as tar sands and bituminuous 
sands, have been found throughout the world, but the largest 
deposit, estimated at 900 billion barrels of pide bitumen, 
lies in the northeastern part of the Province of eect a 
Canada (Baughman, 1978). Despite the fact that this reserve 
has been known for many years, it is only recently that tar 
sands have assumed commercial importance. Difficulties 
encountered in mapping, exploration and processing have 
hampered commercialization of tar sands. 

In the present commercial oil sands plants Such as 
Syncrude and Suncor in Ft. McMurray, Alberta, strip or open 
pit mining has been used. This technique involves the 
removal of overburden and mining of tar sands using either 
bucketwheel excavators or draglines and transportation to 
the extraction plant by a series of long conveyor belts. 
Using the 'hot water' process developed by Dr. Karl A. Clark 
from the Alberta Research Council (Puttagunta et al., 1977), 
the bitumen is then extracted and upgraded. The remaining 
sand is sent to the tailings disposal. Storage of the 
tailings poses a major problem for a surface mine operation. 

The oil sands formations range up to 60 meters in depth 
with overburden depths up to 700 meters (Humphreys et al., 


1977). The oil content of the sand and the thickness of the 
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overburden limit the depth Bi the deposit for which strip 
mining techniques can be used economically. At present, only 
10% of the reserve is considered mineable by surface methods 
(Baughman, 1978). Recovery of the remaining oil sands 
deposits will require either subsurface mining or in-situ 
techniques. 

The principle of in-situ methods is to remove the 
bitumen only, thus eliminating some of the problems 
encountered in the disposal of sand. Of these methods, 
present studies have centred around forward/reverse 
combustion and steam injection. Both methods have undergone 
pilot scale testing in Alberta (Puttagunta et al., 1977), 
but the steam injection method has received the most 
attention and is considered to be exploitable in the near 
future. In this method, high pressure steam is injected into 
a central well. The addition of heat through this method 
into the oil sands formation lowers the bitumen viscosity. 
This enables the bitumen to be pumped from a group of 
producing wells surrounding the central one. 

In ofder sto foptimize *thevextract ion ‘of "Oil; it “rs 
highly desirable to obtain good communication between the 
central and producing wells. One method for accomplishing 
this is to pump a liquid under high pressure into the 
formation in an attempt to produce horizontal fractures. 
However, factors such as the orientation of bedding planes 
and in-situ stresses result in a poor success rate with this 


method. A well is defeated when vertical fractures to 
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formations above or below the deposit are created. Hence a 
controlled method, such as introducing a high pressure water 
jet to produce a series of horizontal holes radiating 
outward from the well bore is of special interest. This 
method could be used to direct the fracturing process in 
establishing a better communication between wells. In so 
doing, the injection wells can be placed further apart 
providing an economic advantage. Other possible applications 
of this method include shaping electrodes for electrical 
heating and increasing the available surface area in the 
bore-hole for better heat transfer. 

The capability of a high speed water jet to fracture 
and penetrate Substances like rock and coal has been amply 
demonstrated (see for example Singh et al., 1974, 
Erdmann-Jesnitzer et al., 1980). In order to investigate the 
feasibility of jet cutting in oil sands, a correlation 
relating the depth of cut to jet parameters is required. 
However, little information pertaining to the depth of 
penetration of water jets in oil sands is available. Some 
earlier laboratory studies and field tests to determine a 
correlation have been carried out by Gilpin and Gates 
(1981). From the laboratory investigation, a correlation 
relating the depth of penetration to jet parameters (nozzle 
diameter and dynamic pressure) and total jetting time was 
developed. The results predicted by this correlation agreed 
closely with those of the field test. However, neither the 


correlation nor the field test provided sufficient 
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information to predict the Gi cdnare depth of the penetration 
or the cutting time. Both of these quantities are important 
for evaluating a jet cutting operation. 

In an attempt to obtain estimates for the ultimate 
depth of penetration, some test results of the 
characteristics of unconfined water jets by Yahiro et al. 
(1980) were incorporated into the first model (Gates and 
Gilpin, 1982). Yahiro et al. measured the average centerline 
pressure of a high pressure water jet as a function of the 
distance from the nozzle ooh various exit conditions. These 
exit conditions were: water jet into air, water jet into 
water, and water jet with a concentric air jet into water. 
Using the criterion that no cutting of oil sands occurs when 
the jet dynamic pressure at the nozzle falls below 1.2 MPa 
for standoff distances greater than 100 nozzle diameters 
(Gilpin et al, 1981), three estimates of ultimate depth were 
obtained. The field data were compared with these 
predictions and it was found that the field data fell below 
the curve for the water jet into air but above that for the 
concentric air water jet into water (Gates et al., 1982). It 
appears then, that the maximum depth probably lies between 
the curve for a water jet into air as an upper bound and 
that for a concentric air-water jet into water as a lower 
bound. A comparison between these estimates based on the 
lower bound and two cutting theories (Crow (1973), Hashish 
and duPlessis (1978)) was also carried out (Gates et al., 


1982). The results of the comparison is demonstrated in 
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figure 1.1. From this figure, it can be seen that Crow's 
theory substantially underpredicts the maximum depth of 
penetration and also shows a linear dependence upon pressure 
while the actual case is nonlinear. A better agreement with 
the experimental data was found with the theory of Hashish 
and duPlessis by adjusting a parameter Cf which is a 
hydrodynamic coefficient of friction. However, the cutting 
theories were developed for a traversing jet rather than a 
Stationary one and extension to the stationary jet case is 


questionable. 


1.2 Present Investigation 

As can be seen from the above discussion there is no 
theoretical method or sufficient experimental data from 
which a prediction of the ultimate depth of penetration can 
be made. It is desirable then to have results similar to 
those of Yahiro et al. for a confined jet as this would 
provide a more realistic model of the actual situation than 
an unconfined jet. To this end an experimental apparatus was 
built to simulate a confined jet in a down-hole jetting 
geometry. The flow of the jet was investigated qualitatively 
by flow visualization and quantitatively by measuring the 
pressure of the jet on a simulated cutting face. From this 
experimental data an empirical model to predict maximum 
depth of penetration and cutting rates was developed. The 
predictions of this model were then compared to some limited 


field test data and to a correlation based on Gilpin (1981). 
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2. Experimental Apparatus and Procedures 


2.1 Apparatus 

A laboratory model was constructed to simulate a 
down-hole jet cutting geometry like that of figure 2.1. The 
experimental apparatus is shown schematically in figure 2.2 
and a photograph of the apparatus is provided in plate 2.1. 
The model consists of a vertical and an horizontal section 
Supported in a steel tubing frame. The test model was 
constructed with plexiglass in order to study the flow 
visually. 

The vertical section consists of a 0.3 meter square 
tube, 2.5 meters high which is open at the top and closed at 
the bottom. Drainage of the spent water at the bottom of 
this vertical section was controlled by two valves. A square 
window of 0.2 meter was constructed on one side of the 
vertical section to provide access to the nozzle assembly. 

The horizontal section representing the hole pierced by 
the water jet was modelled by a circular tube. Because the 
Size of the actual hole produced for given jetting 
conditions is not known several sizes of tubing were used 
ranging from 38 mm to 102 mm in diameter. A moveable piston 
carrying a load cell slides along the length of the tubing 
simulating the depth of the hole cut. The piston is 
connected to a long tubular rod which is held against the 


frame using a tapered clamping device. Hence, the piston can 
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piston 


nozzle assembly 


Plate 2.1 A photograph of the experimental apparatus used 
to model a down-hole jet piercing geometry. 
In this photograph a 38 mm diameter horizontal tube 
is in place and the piston to nozzle distance 
is 0:9 meter. 
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be fixed at any position along the tube. 

In addition to the solid walled tubes, some experiments 
were carried out uSing a porous walled tube. The motivation 
for these tests came from field test observations (Gilpin 
and Gates, 1981) in which it was noted that spent water 
returned through previously drilled holes and cracks in the 
oil sands formation. Consequently, the assumption of a solid 
wall appears to be too severe. To simulate the porosity, a 
large number of holes were drilled in a 51 mm diameter tube. 
Five percent of the total surface area of the tube was 
drilled out uniformly over the entire length of the tube, 
with holes of 6.35 mm in diameter. A 0.2 meter square box 
was built to enclose the porous tube to prevent splashing of 
the water when the porous test was done. Spacers supported 
the porous tube inside this box. This box was open at the 
back allowing the water to drain away. This arrangement can 
be seen in plate 2.2. 

A high pressure three cylinder positive displacement 
pump driven by a 30 kW electric motor generated a 
pressurized flow of water which was conducted through a 
flexible hose to the nozzle assembly. The volume flow rate 
of the water was determined from the rotational speed of the 
pump. Pump speed was measured utilizing a magnetic sensor 
and a magnet attached to the hub of the pump drive pulley. 
The resulting pump calibration curve is shown in figure A.1 


in the appendix. 
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Plate 2.2 A photograph of the apparatus arranged for 


porous wall tube tests. The square horizontal 
box enclosing the porous tube collects the 
water escaping through the pores. 
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2.2 Nozzle Assembly 

The nozzle assembly mounted in the vertical section 
opposite the horizontal tube was designed to absorb the high 
thrust force of the jet and to allow for fine adjustment of 
the nozzle position and direction. This assembly, using 
several set screws and washers, was capable of adjusting the 
jet in three degrees of freedom: yaw angle and vertical and 
horizontal displacements. A photograph of this assembly can 


be seen in plate 2.3. 


2.3 Nozzle Geometry 

The nozzles used in this experiment were made of 
Stainless steel with a tungsten carbide insert. As no 
information pertaining to the geometry of the nozzle was 
available from the manufacturer (Hydro Manufacturer Inc.), 
the internal geometry of the nozzle was determined by taking 
a mold with ERTV molding rubber. This molding has a volume 
shrinkage of less than one percent. The nozzle profile was 
then determined by placing the hardened mold ina 
shadowgraph and recording the profile coordinates. The 


nozzle profile is shown in figure 2.3. 


2.4 Instrumentation 


2.4.1 Load Cell 


To measure the pressure that the water jet exerts ona 


material surface, a load cell was constructed and mounted in 
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Plate 2.3 A photograph of the nozzle Support and 
alignment assembly. 
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Figure <2'.3 


D: Nozzle diameter 


A schematic diagram of the nozzle geometry. 


15 


i 
2 
» 4 
eae 
{ uy 
i ' 
I 
iy 
| i 
1 
\ 
7 
f d : 7 a 2 an 
c - ii ' 
, 
‘ Pe 1) 7 On 
ae 3] 
. a” , ‘ 
- a a y as es '? \ 
ie usemasp eit 
a | : : 


: 7 nal re 
- ip 


16 


the piston surface. All the pistons were made to accommodate 
the same load cell. The piston and load cell arrangement can 
be seen in plate 2.4. 

The load cell is a 25.4 mm diameter by 2.54 mm thick 
steel diaphragm with a 6.4 mm rim as shown in figure 2.4. 
Foil strain gauges were mounted on the side of the diaphragm 
that was not exposed to the water. Four gauges were used to 
give a full bridge configuration. The arrangement of the 
gauges, Similar to Werner (1953) as shown in figure 2.4, 
provided temperature compensation. With a uniform pressure 
applied to the diaphragm, the middle gauges (2 and 3) are 
subjected to tension while the outer gauges (1 and 4) sense 
compression. The strain measured from these gauges is 
directly proportional to a uniform pressure applied to the 
diaphragm. An experimentallly measured value of the first 
resonant frequency of the plate showed a value of 20 kHz. 
This was done by dropping a small ball bearing onto the 
middle of the load cell and observing the trace of the 
strain gauge bridge output on an oscilloscope. 

The load cell was statically calibrated at room 
temperature with compressed air as shown in plate 2.5. The 
relationship between the applied pressure to the measured 


electrical output is shown in figure A.2 in the appendix. 


2.4.2 Signal Analysis 
A schematic diagram of the load cell signal analysis 


procedure is shown in figure 2.5. The signal obtained from 
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Plate 2.4 


load cell 


A photograph of the four pistons used in these 
experiments. Also shown is the load cell about 
to be mounted in the 102 mm piston. 
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Figure 2.4 A schematic diagram showing the load cell 
geometry and the positioning of the strain gauges. 
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Plate 2.5 A photograph of the apparatus used to 
calibrate the load»cell. 
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the load cell after amplification had a very large 
fluctuating component and a small steady offset. A sample 
trace of the signal from the oscilloscope is shown in plate 
2.6. Because of the nature of the pressure signal, it seemed 
more appropriate to examine the signal in more detail than 
only recording peak and average values. The discrete nature 
of the signal suggested that the impact of the water jet 
could be viewee as a stochastic process. The pressure signal 
is therefore analyzed to provide a frequency of occurrence 
distribution. This was accomplished by recording the signal 
with an HP 3721A signal correlator to provide the frequency 
and cumulative frequency distributions. From frequency 
distribution, the probability density function is determined 
by normalizing the frequency distribution to produce an area 
under the curve of one. The area underneath the probability 
density function is known as the cumulative density 
function, and gives the probability that a signal is below 
or above a certain amplitude. 

The pressure signal was sampled at a rate of 3 kHz (the 
fastest available), and each sample is placed into a bin 
according to its magnitude. In probability density 
measurements the sampling sequence stops when the number of 
samples in any one bin reaches a pre-set number (in the 
present case, 32x1024). In the probability integral 
measurements, processing ceases after about 4 times this 


number of samples have been taken. 
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Plate 2.6 An example oscilloscope trace of the load 
cell signal. 


Vertical scale: 50 kPa/cm 
Horizontal scale: 5 millisec/cm 
Nozzle diameter: 1.6 mm 

Nozzle pressure: 10 MPa 
Standoff distance: 0.5 meter 
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In addition to the frequency distribution the average 
value of the signal was also measured. This was done using a 
voltage to frequency converter and a counter. The input 
pressure signal was converted to frequency and the average 
frequency was then obtained over a fixed time interval (10 


seconds). 


2.5 Test Procedure 

At the beginning of each test the tubing was levelled 
and the jet was centered by adjusting the set screws in the 
jet assembly so that the jet hit at the center of the face 
of the load cell. It was then assumed that the jet hit the 
center of the face of the piston along the length of the 
tube. The jet offset due to gravity was estimated to be 
approximately 0.2 mm for a standoff distance of 1 meter. 

Before any readings were recorded, the load cell was 
checked for seating by running the pump on and off several 
times. In doing this, the load cell would be properly seated 
and the bridge could be balanced. Subsequently, any strain 
measured would be due to the pressure of the water jet. The 
water temperature for all the experiments was kept constant 
atei2 2G 

During the test, for a given nozzle diameter, tube 
diameter and flow rate, the average pressures for each 
Standoff distance were obtained. Plots of probability and 
cumulative density functions from the correlator were also 


recorded. To prevent damage to the load cell, the minimum 
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standoff distance was kept a6 om em. 8Eheitoad cell was 
recalibrated regularly, particularly after running the 
smallest nozzle diameter because of its high pressure, and 
whenever a new piston was used. This was done to ensure no 
permanent deformation of the diaphragm resulted from the 
previous run and the same calibration curve could still be 


used. 
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3. Presentation of Results 


For a given nozzle diameter, dynamic pressure and tube 
Size, plots of probability and cumulative frequency 
distributions were recorded for each nozzle to piston 
(standoff) distance. In addition average pressures were 


recorded. The range of variables used were: 


A. Solid wall 
Nozzle sizes: 3.2 mm, 2.4 mm, 1.6 mm. 
hiubesdiavest £O2 mm. >4/6enmy Siimm, 38 mm. 
Pressures : up to 35 MPa. 


Standostsdierances:<033um todl.6é m. 


B. Porous wall 
Nozzle sizes: 3.2 mm, 2.4 mm, 1.6 mm. 
Tube wdie.\ 7.5) mm. 
Porosity: }9%. 
Pressures : up to 35 MPa. 


Standopt drstance 20.3imvtonds 7ceme 


In addition to making quantitative measurements of the 
pressure the water jets were also studied qualitatively by 
flow visualization. Since the flow observations are very 


useful in interpreting the quantitative data, they shall be 


presented first. 


25 


yi mn is a in 


cs 
20 0 ee ar 


: _ Ke ry r 


Me ea ia oe, ue 


ae dl 


a 
st a He 4 


ae ae 
. 192390 Ga- Siasen nevi g 


ae 
7 


i(ibhes @l ,@oe 


Sabdoe. saa nooks 


26 


3.1 Flow Visualization 
The flow visualization was accomplished by high speed 
photography. A high speed movie camera (HYCAM) was used to 
record the flow of the unconfined and confined jets. The 
camera, running at 5000 frames per second, was positioned 
perpendicular to the jet center-line at a point 0.5 meters 
from the nozzle. Illumination was provided by tungsten 
lamps. From these flow observations it could be seen that 
several distinct flow regimes could be defined. These are: 
1. Unconfined jet. 
2. Confined jet with solid wall: 
a. no slug flow regime 
b. intermittent slug flow 
Cc. permanent slug flow. 
3. Confined jet with porous wall. 
The features of each of these flow regimes will now be 
described. 
1.  Unconfined jet. 

Before examining the more complex features of the 
confined jet flow, the flow characteristics of an 
unconfined jet were investigated. A sequence of frames 
from the high speed movie for an unconfined jet is shown 
mo plate 3.1. The direction of the water jet 1s from 
right to left. As can be seen from this series of 
frames, the jet has a definite structure that is 
repeatable.) The basic structure 1s that of a “slug” jor 


"packet" of water trailed by a large number of small 
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Direction of flow is 
from right to left 


Nozzle diameter: 2.4 mm 

Nozzle pressure: 2 MPa 

Standoff distance: 0.5 meter 

Width of frame: 23cm 

Time between frames: 2 milliseconds 


Plate 3.1 A sequence of high speed camera frames showing 
the flow of an unconfined continuous jet. 
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droplets. The reason for this structure and its 
repeatibility can be traced to the manner in which the 
jet is produced. 

Recall that the jet 1S produced by a positive 
displacement pump having three pistons. On each 
discharge stroke as the piston moves from bottom dead 
center to top dead center its speed varies from zero to 
some maximum value and then back to zero. Consequently a 
Similar speed variation occurs in the water jet as it 
leaves the nozzle. The high speed water produced near 
mid-stroke very rapidly overtakes the slower moving 
portion of the jet produced at the beginning of the 
stroke forming a "slug" or."packet" of water. The slower 
moving portion of the jet produced near the end of the 
jet trails even further behind the slug in a series of 
small droplets. 

This modulation of the water jet was clearly shown 
by a spectral analysis of the pressure signal and 
correlated exactly with the pump rotation rate. 
Confined jet with solid wall. 

Observations of the flow when the jet was confined 
within the solid wall tubes suggested that the flow be 
divided into three regimes. These are: 
a. No slug flow regime. 
In this flow regime no interference occurs 
between the incoming jet and the spent water. The 


water jet hits the piston and flows radially across 
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the piston face towards the tube wall. As it hits 
the tube wall, the water jet turns 90 degrees and 
flows along the tube wall towards the incoming jet. 
At a small standoff distance, spent water returns 
along the tube wall concentrically to the incoming 
jet. As the standoff distance increases, the spent 
water from the upper part of the tube falls along 
the tube wall and joins the flow of spent water from 
the lower part of the tube. A sketch of the latter 
type of flow is provided in figure 3.1. This flow 
regime was found in the 76 mm and 102 mm tubes 
almost over the entire range of standoff distances. 
Intermittent slug flow. 

As the standoff distance is increased some 
interference between the spent water and the 
incoming jet occurs. This interference can be caused 
by either rapid spreading of the jet so that it 
contacts the tube walls or by large flow rate so 
that the spent water can not drain fast enough. The 
result of this interference was to cause a "hold up" 
of some of the spent water so that a slug of water 
forms in front of the piston. The presence of the 
slug at this stage is not continuous, but depends 
upon which portion of the jet is impinging on the 
piston. A sketch of the flow can be seen in figure 
Sees 


Permanent slug flow. 
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As the standoff distance is increased further, 
the accumulation of the spent water from the 
intermittent slug flow increases to the point where 
the presence of the slug becomes continuous. At this 
point the incoming jet no longer strikes the piston 
Surface directly, but instead is dissipated at the 
front or upstream face of the slug. Further 
increases in the standoff distance results ina 
lengthening of the slug with the upstream face 
remaining approximately in the same position. A 
sketch of the flow in this case is provided in 
LiGune (3.3): 

Confined jet with porous tube. 

The holes in the porous tube allow the spent water 
to drain from the tube preventing any accumulation of 
water at the piston face. The spent water flows radially 
from the piston towards the tube wall. Most of the water 
flowed directly out of the tube through the pores when 
it reached the tube wall. However, some of the spent 
water changed direction as it hit the wall and flowed 
towards the incoming jet. The pores permitted this spent 
water to drain from the bottom of the tube without ever 
reaching the vertical section. A sketch of this type of 


flow can be seen in figure 3.4. 
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Figure 3.1 


A sketch illustrating the flow within the 
tube when no interference occurs between 
the incoming jet and the returning water. 
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Figure 3.2 A sketch illustrating the flow within the 


tube when interference between the incoming 
jet and the spent water begins to occur. 


air water 


Figure 3.3 A sketch illustrating the flow within the 
tube when the presence of the slug has 
become continuous. 
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Figure 3.4 A sketch illustrating the flow within the 
tube of a porous wall. 
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3.2 Quantitative Results 


3.2.1 Frequency distribution 

Plots of the probability and cumulative density 
functions for various distances and sets of parameters were 
recorded on a X-Y recorder and later digitized. Some 
examples of these results are provided in figure 3.5 to 
3.11. The results are recorded in Gates and Toogood (1983). 

From figure 3.5 it can be seen that for small standoff 
distance the mean and variance of the probability density 
function are large. As the distance from the nozzle 
increases, the probability density function becomes narrower 
indicating a decreasing variance. The mean of the pressure 
Signal also decreases as the standoff distance increases. 
This general trend is observed for smaller tube diameters as 
Shown in, figures: 3.1.60 3.9. The shape of the probability 
density function for large tube diameters rises and drops 
Oftmouaduallys (figures 3.5 to 3.7). However, for smaller 
tube diameters where the slug forms, the shape of the 
probability density function behaves differently as 
indicated in figures 3.8 to 3.9. The probability density 
function rises and drops off immediately without any gradual 
tailing off as observed in the non-slug case. The variance 
of Ehrs slug case, 1Sssmall. This TSsedue to the tact that the 
slug of water acts as a cushion in front of the piston hence 
dampening any pressure fluctuations. For a higher dynamic 


pressure, the general trend of the probability density 
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Solid wall (102 mm) 
Go¥ =!) 3% 22mm 

Pg area 3 MPa 
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Pressure at the cutting face (kPa) 


The probability and cumulative density 
function *of ‘the “load cell signal for 
increasing values of the standoff distance 
(D=102) mm, Pi =23 MPa)’: 
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Pressure at the cutting face (kPa) 


Figure 3.6 


The probability and cumulative density 
function of the load cell signal for 


increasing values of the standoff distance 


(D=102 mm, Pp, =20.5 MPa). 
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Solid wall (76 mm) 
do =s 3% 2am 
=) 1%. 3.MPa 
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Pressure at the cutting face (kPa) 


The probability and cumulative density 
finctaoneof. the WVoads cell signal for 
increasing values of the standoff distance 
(D=76 mm, P,=1.8 MPa). 
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Figure 3.8 


Solid wall (51 mm) 
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Pressure at the cutting face (kPa) 


The probability and cumulative density 
function of the load cell signal for 
increasing values of the standoff distance 
(D=5) mm sold, P,=2.2 MPa). 
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Solid Wall (38 mm ) 
3.2 mm 
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Pressure at the cutting face (kPa) 


The probability and cumulative density 
function of the load cell signal for 
increasing values of the standoff distance 


(D=38mm, Py=2.1 MPa). 
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Pressure at the cutting face (kPa) 


The probability and cumulative density 
function of the load cell signal for 
increasing values of the standoff distance 


(D=38 mm, P=15.6 MPa). 
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Porous wall (51 mm) 
do = 1.6 mm 
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Pressure at the cutting face (kPa) 


The probability and cumulative density 
funet von of themtoad cell signal tor 
increasing values of the standoff distance 


(D=5.1 mm porous, P,=21.4 MPa). 


‘aa a 
Vite wt eee 
ay * 
| it 
A il - Zi * rr > 7 
: , 
‘s , na > A a a 
~ 7 14 7 
> wae jl 
‘ : cn 
wd q 
1 , 
=] ig ; 
’ : tn Aa 
4 “1 : 
; 7 
: 2 3 ae See ec 
al ' : bs e : 5 i 7 i Se ga ae ;, = 
= et ee ™™, “Sa a; 
sete oil "7 = - ih _o ' +> sh ew Ye 
iy," ; . be) . ; 
~* ee oy a : L oe, 
_ ™ p 
- 4 |* 
; i 
4) 


” % i a 
/ a Vi, - | 
Me ¢ 
, i” eli qe Sy ba 
. ne 3 
eae ye pas ca a) m4 
ra x Te 
, oh ¢G , i 8 y 
p « " : 
| : : vi i - ° a “¥ 
es ats. a 
t u é rt 
‘6 7 - 
HT = Be 
e 7 
wer 
7 f fay ; 
7 ff 
' a Fay 
a : uw ry 
mt ; s i ~ 
i. hoe fee 
‘ i = : if i " 
2 y ; 7 
a 7 
J ¥ 
{ hid ae 2 7 
ete, ' ak S i 
— } ¥ & _ 
* iS Lie P Pa . “ 
ms | aloe 
i a a i Tih ee : p, a 
OM ie. 
a ae . ] cy 
> J 7 gow A 
: se @ Hel i. £ 
. a) ie “ , 
. 4 _ 
t bale A ihe Ven’: oA yoy ae mM 


. : oe 
bert) se e ree 9a Por 
tp ms 4 wi ais 


“ 
7 my 


: ay. 
5 ae L2rob ry vi ins 7 sf : f / : 
“, . 712 PEnSi se. . 5LRe 5 i 

7 "ong ith ‘eokn ado ‘gt? 1% 3 


e era) 3 


4 
7 


40 


FUNCTION sti itenolds. That 16 asi the distance from the 
nozzle increases, the variance and mean of the pressure 
Signal decrease. The only difference is the magnitude of the 
mean and the variance. For higher dynamic pressure, the mean 
and the variance are much higher than those of for lower 
dynamic pressure as indicated in figures 3.6 and 3.10. 

The same general trend mentioned earlier is also 
observed for the probability density function for the porous 
wall tube as shown in figure 3.11. From the flow 
visualization, no slug formation is observed due to the 
pores of the wall. This is verified from the shape of the 
probability density function where all of the curves drop 


OLr Qradualtly. 


3.2.2 Average pressure 

The average pressure waS non-dimensionalized using the 
dynamic pressure and plotted against the standoff distance 
to nozzle diameter ratio for each tube diameter. The results 
are summarized instne Log-log plots, in f1gures. 3.12 to 3.46. 
Most of the data falls onto a straight line. At high 2/d> 
the curve levels off for the 38 mm and 51 mm solid tubes. 
From flow visualization, the distance where the graph starts 
to level off indicates the slug formation. The measured 
pressure is independent of the standoff distance once the 
slug. forms. Thus verdne Othe wact nat lie distance where 
the slug forms stays at the same position and all of the jet 


energy is dissipating at that distance. Further increase in 
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Figure 3.12 A plot of the non-dimensionalized average 
pressure against the non-dimensionalized standoff 
distance for the 102 mm diameter solid wall tube. 
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Figure 3.13 A plot of the non-dimensionalized average 
pressure against the non-dimensionalized standoff 
distance for the 76 mm diameter solid wall tube. 
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wall (D=51 mm) 
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Figure 3.14 A plot of the non-dimensionalized average 
pressure against the non-dimensionalized standoff 
distance for the 51 mm diameter solid wall tube. 
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Porous wall (D=51 mm) 
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Figure 3.15 A plot of the non-dimensionalized average 
pressure against the non-dimensionalized standoff 
distance for the 51 mm diameter porous wall tube. 
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Solid wall (D=38 mm) 
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Figure 3.16 A plot of the non-dimensionalized average 
pressure against the non-dimensionalized standoff 
distance for the 38 mm diameter solid wall tube. 
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the standoff distance results in a lengthening of the slug 
with the upstream face remaining approximately in the same 
position. Also from these figures, it can be seen that the 
dependence of P ave/Py with Z/d»o is independent of the tube 
Size provided no slug forms. This is demonstrated by 
Superimposing figures 3.12 to 3.16. All the data fall onto 
the same straight line and start to deviate only when the 
slug forms. | 

The results of the average pressure for the porous wall 
tests are plotted in figure 3.15. No slug formation was 
observed. This is confirmed from figure 3.15 where all the 
data falls onto a straight line, with no levelling off of 


the curve as in the solid wall cases. 
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4. Analysis of Results 


Flow visualization observations and pressure signal 
records revealed that the effect of the water jet on the 
piston surface consists of a series of discrete impulses 
rather than a continuous but unsteady force associated with 
a coherent jet. It 1S suggested here then that an 
appropriate approach is to consider the jet impact to be a 
stochastic process and to develop a probabilistic model for 
predictions of ultimate depth of penetration and the cutting 
time required to reach that depth. In the model proposed 
here, it is first assumed that there is some critical 
pressure (Pc) that must be exceeded to cause failure at an 
oil sands surface. For each set of values of the standoff 
distance, nozzle diameter, nozzle pressure and tube 
diameter, the measured cumulative density function will then 
provide the probability thatetheveriticalspressureryrs 
exceeded, ie. 

@ (p>Pc) 
where @: probability that the pressure exceeds Pc; 
Pe: critical pressure; 
Dp: (pressure at the cutting face. 
Holding all parameters constant except standoff distance, a 
plot of this probability against standoff distance can be 
made. Now if it is further assumed that @(p>Pc) must exceed 


some value @* for a non-neglible cutting rate to occur, the 


47 


5 
¥ ‘ © ‘ ¢ - 
oe I u e ox “ 
ea - 
, a : 
= A 
- y 


a 
sje2 fe ete “uid Bots sxe abet 308 7 
nn . a fie sok 7 ail + oF 
ty yas Eo etree some 78" oT | 
hi hy 5 a3 ah See Svouns.angds ‘Si aade vada ‘ 
. end —eeteegpis af rl \aet Preet 
: : : ope Ser t- a” 


al yi . ae) 


er at hianbe o4 Bf Ade tis 916s i903 ag 


Piper epi idedoy¢ « Golaveatl Gf Das ea s3074, Jiaiaed 


a oF Pee 7 Tsao SIRT iL. 36 eno fast 
ue 


tuseeota twleh ots nt gees sete wee oF bsatupih 


ai : f a i - in i! 
a oe nig = y y 
aL | eels em Ome b Saath 68 4 2 | ae sf it 


ulisl. sepst ty Sebdeoy ‘ae aa inky a) aayaeotg 


iad 
y 7 ; i? 
y Meee eyes PRO > sandaines eisaba Bip aie 


| ror hese a Ay .) “y,) 

Liga nord nee waeah av (Pola! Latenng ‘sanwaton naa \redsbeit 
- 

et. wy sce, Levtidian see amis gat ents ont ab. a&ivOI 1¢ 


% Oe) de iY tw La fon 2) 


it * 5 | {3 7 7 & SV Wa * a Ba Aaee ) ‘Ys ix. a5 230 co “a tte 


oe 


fi3 ie het pants su 


i.e smnsuets iebnsta sgpune A ais 
eG Nie. . 
40 sik ak sauntei® hinge ‘38 


Ai =_ 


tonne: — Daisice bet 


+ edi - £ ney 315t 


a 


e iy 


48 


ultimate depth of penetration can be determined. In 
addition, if a relation between the penetration rate and the 
probability is known, an estimate of the time required to 
reach the maximum depth can be obtained. The details of this 


model will now be presented. 


4.1 Prediction of Maximum Depth of Penetration 

First, an estimate for a critical pressure is required. 
This estimate may be obtained by examining the experimental 
results of Chau (1980). Chau observed that no penetration of 
oil sands occurs at standoff distances greater than 200 mm 
when the dynamic pressure of water jet falls below 1.2 MPa 
from a 1.6 mm diameter nozzle. Although the dynamic pressure 
at the nozzle outlet is known, no information on the 
pressure at the cutting face was obtained by Chau. To obtain 
an estimate of this, the above conditions were duplicated 
including the use of the same nozzle. Plots of the 
probability and cumulative density functions were recorded 
and are shown in figure 4.1. Using the assumption that no 
significant cutting of oil sands occurs when the probability 
{hor the jreceure exceeds the critical value falls below 
5%, a critical pressure of 70 kPa was determined (see figure 
aerley, 

With this value for a critwveal pressure, ithe 
probability that the pressure at the face of the piston is 
greater than Pc for various distances was obtained from the 


measured cumulative density functions. These probabilities 
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@(p>Pc)=0.05 


Py=1.2 MPa 
do=1.6 mm 
Z=200 mm 


Pressure at the cutting face (kPa) 


Figure 4.1 Probability and cumulative density function 
using the same nozzle as Chau (1980). 
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were obtained based on non-Slug experimental data. It was 
assumed that the depth of penetration reached its maximum 
depth when the slug formed. These probabilities are provided 
in figures 4292 tor Ge Toe"each™set ot data’ a curve of the 


form: 


ef = exp(-A’(Z/do)*) (1) 


was fitted where 

Z : depth of penetration; 

de’ 2) nozzievdirameten: 
and \} is a parameter that depends on the dynamic pressure 
P,, nozzle diameter do, and tube diameter D. 

The curve of this form was chosen because the rapidly 
decreasing probability with distance from the nozzle appears 
to be exponential in nature. The quadratic form of Z/do was 
chosen instead of the linea one because it would be expected 
that @(p>Pc) would be near one for some distance before 
beginning to decrease. 

The values of > in equation (1) for the various test 
conditions were determined by the method of least Squares 
and the results»ares-plottdd-in-fa.gurer4.7..-From@thiseiigure, 
it can be seen that for the range of the current tests, A 
depends most strongly upon the dynamic pressure at the 
nozzle. This parameter A is observed to be approximately 


inversely proportional to dynamic pressure. A relationship 
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Figure 4.2 
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piston face exceeds the critical value for 


some combinations of nozzle diameter and 
pressure in the 102 mm diameter tube. 
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b> 
Solid wall (D=76 mm) 

™M Pa= 1.8 MPa do=3.2 mm A=0.485x10-4 

A P4g= 2.6 MPa do=3.2 mm \=0.333x10-4 
4 & Pa=11.4 MPa do=1.6 mm X=0.109x10-! 
“| *#P,=19.7 MPa do=1.6 mm X\=0.616x10-3 

Z P,=34.6 MPa do=1.6 mm X=0.306x10-5 

€(p>Pc) =e 78! 2/de) * 


Figure 4.3 
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The probability that the pressure at the 
piston face exceeds the critical value for 
some combinations of nozzle diameter and 
pressure in the 76 mm diameter tube. 
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Figure 4.4 
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1000 


piston face exceeds the critical value for 


some combinations of nozzle diameter and 


pressure in the 51 mm diameter solid tube. 


Probability (Pressure > P critical) 


Porous wall (D=51 mm ) 


Figure 4.5 


do=3.2 mm A=0.322x10-4 
do=3.2 mm \=0.264x10-4 
do=1.6 mm A=0.607x10-5 
do=1.6 mm A=0.281x10-5 
Go=1.6 mm A=0.190x10-5 


200 400 600 800 1000 


Distance/Nozzle diameter 


The probability that the pressure at the 
piston face exceeds the critical value for 
some combinations of nozzle diameter and 
pressure in the 51 mm diameter porous tube. 
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Probability (Pressure > P critical) 


Solid wall (D=38 mm) 


MPa do=3.2 mm »\=0.813x10-4 
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@(p>Pc)=e ~*'2/do)* 


Distance/Nozzle diameter 


The probability that the pressure at. the 
piston face exceeds the critical value for 
some combinations of nozzle diameter and 
pressure in the 38 mm diameter tube. 
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Tube Diameter (mm) 


oOO1-6 
Ov 16.2 
4. 50:8 
+ 38.1 
x 50.8 (porous) 


| 32mm Nozzle 


ie 


2.4mm Nozzle 


Se 


a7 sO 1.6mm Nozzle 


Jet dynamic pressure (MPa) 


Figure 4.7 The variation of the parameter A with 
dynamic pressures for different tube and 
nozzle diameters. 
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between A and dynamic pressure was assumed as follows: 


i aye, © 2) 


where a and k were determined by least squares method to be: 


a OT To tx eS 


kK 


1.0986 where Py iS given in MPa. 
Therefore, equation (1) can be rewritten as: 


{raha} 


fF =e (3) 


It is further assumed that no cutting occurs when the 
probability that the pressure exceeds critical pressure 
falls below 5%. Substituting @ orobability of SY into 
equation (3), the resulting prediction of maximum depth of 


penetration becomes: 


Ze = ¥(-(do?Pytin C)/ay (is) 


where © =s0%705 


This equation can be simplified approximately EO: 


Zult = l65*dadoyvel2i (5) 
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where Zult is in mm; 
ous tim mm: 
Pp fis tim MPa. 

The dynamic pressure of the water jet is related to the 


volume flow rate by: 


Moms ce : 
2 3 x 10 ( MPa ) (6) 


where p : density of water; 
Q : volume flow rate of wates jet. 


Therefore, equation (5) can be also approximated by: 


PSO ee 
Zult = 1 RAO do (7) 


where Ovis an m?/sec* 

avis * in mins 

Zult is in meter. 
The above relation shows the dependence of the ultimate 
depth to the flow rate and nozzle diameter, and is validyvion 


conditions where no slug forms. 


4.2 Prediction of Cutting Time 
The time to reach the ultimate depth is another 
parameter of interest for an operator of this water jet 


instrument. It is required to predict the water requirements 
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for this (type OLsapplication. 

Numerous cutting relations might be proposed, however, 
one assumption is of direct proportionality between the rate 
of cutting and the probability that the critical pressure is 
exceeded. A better cutting relationship could be 
investigated; however, before any more complex treatment of 
this model is carried out, a simple model approach will be 
done. Hence the rate of penetration based on the direct 


proportionality takes the form: 


€@Z/dt =.Be (8) 


where dZ/dt : cutting rate, and 

B : a parameter which has dimensions of length per 
unit,.time. 
This parameter B can be interpreted as the initial cutting 
rate, because at time t=0, the probability of the pressure 
exceeding critical pressure Pc (provided P,>Pc) is equal to 
one. Substituting equation (3) into equation (8), and 


integrating yields 


(9) 


which gives the time necessary to reach a given depth of cut 


Z. In particular, with the ultimate depth given by equation 
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(4), the time required is: 


2 
| me (10) 


The value of the parameter B was determined from the 
laboratory data of Chau (1980). Only the data obtained from 
1.6 mm nozzle diameter could be used, because the dynamic 
pressures for other nozzle diameters were rather low. The 
value of B in general would be a function of standoff 
distance, dynamic pressure, nozzle diameter, and material 
properties. The method to obtain this parameter B is 
jllustrated in the appendix  B. The resulting initial cutting 
rates for different dynamic pressures are plotted in figure 
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Figure 4.8 Estimates of the parameter B based upon 


best fits of equation (9) to the: laboratory 
datasot Gnaw 980). 
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5. Discussion of Results 


5.1 Comparisons of prediction and other correlations 

The depth of penetrations predicted using equation (9) 
for 3 different dynamic pressures are plotted and compared 
with some field test data (Gilpin, 1981) in figure 5.1. As 
can be seen from this figure, the prediction of the depth of 
penetration agrees fairly well with the field results. A 
maximum of 8% deviation for the dynamic pressure of 17 MPa 
was observed between the prediction and the field test 
results. For the 12 MPa case the comparison is not so 
promising. However the large scatter in the field data makes 
ameaningful companison @ifiicult,. The scatter in the field 
test results is believed to be due to the orientation of the 
bedding planes. This orientation is not taken into 
consideration in the laboratory correlation developed here. 

A comparison between field test results and predicted 
ultimate depth of penetration for various dynamic pressures 
is shown in f#aqure 5.52. Most cf the field test data are 
below the predicted ultimate depth which is as expected, for 
the ultimate depths were not attained. All of the data 
available from the field tests is for very short jetting 
times (60 seconds). The cutting times predicted for the 
ultimate depth are much longer(eg. 300 sec vs 60 sec). 
Hence, no conclusive comparison can be made on the predicted 


ultimate depth or required cutting time. 
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Figure 5.1 A comparison between depths of penetration 


predicted by equation (9) and the field data 
of Chav (1980) for a 6.4 tim diameter nozzle, 
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Also plotted in figure 5.2 are the predicted depths of 
penetration for a continuous jet using the correlation 
developed for recompacted oil sands (Chau, 1980). The 
results show that the depth of penetration based on the 
probabilistic model is lower than that of the correlation 
developed by Chau. This is not too surprising as first 
Chau's correlation is based on laboratory data for very 
limited jetting times (4 seconds). Second, Chau's 
correlation predicts no limiting depth but states the depth 
should continue to increase indefinitely in proportion to 
the 0.44 power of time. 

The ultimate depth of penetration predicted by Gates et 
al. (1982) is compared with the present correlation for 
various dynamic pressures in figure 5.3. This plot shows 
that the ultimate depth predicted based on a confined jet is 
slightly below that of unconfined jet for dynamic pressure 
below 24 MPa. As the dynamic pressure increases, the 
predicted depth based on confined jet is higher than that of 
unconfined jet. 

The simplified form of the ultimate depth given by 
equation (5) shows the ultimate depth is proportional to the 
nozzle diameter and square root of the dynamic pressure. 
Thus increasing the nozzle diameter or dynamic pressure will 
cause the maximum depth of penetration to increase as well. 
Increasing the volume flow rate will also increase the 


ultimate depth of penetration as seen from equation (7). 
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Figure 5.3 A comparison of ultimate depths predicted by 
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5.2 Influence of some assumptions in the empirical model 

It has been assumed that when the probability that the 
pressure exceeds a certain critical pressure falls below a 
given value, an ultimate depth has been reached. Hence, this 
approach depends on the critical probability chosen. Figure 
5.4 shows the dependence of critical pressure obtained on 
the probability chosen. If the probability that the pressure 
exceeds the critical pressure increases, the corresponding 
critical pressure decreases and vice versa. For example, if 
a critical probability of 30% was chosen, the corresponding 
critical pressure that must be exceeded is found to be 40 
mPaeas shown an figure’ 5.4.88ased onithis critical pressure, 
new experimental constants for equation (3) can be obtained. 
The ultimate depth of penetration predicted using a critical 
probability of 30% was about 25% higher than the one 
predietedvusing al critical probabiiityrvol 5%.) The critical 
pressure depends on the material properties and other 
factors such as the overburden pressure and the orientation 
of the bedding planes in the formations. More data would be 
very beneficial in providing a better estimate of this 
quantity. 

The simple linear relationship assumed here for the 
penetration rate seems to provide a reasonable estimate of 
the time of penetration. Physically, however it is 
unsatisfying since it means that the rate of penetration is 
independent of the pressure at the cutting face. A more 


reasonable relationship for the cutting rate can be proposed 
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by including the effect of the amount by which the expected 
value of pressure exceeds the critical cutting pressure. 

The correlation developed was based on the non-slug 
experimental data. When the slug forms as in the case of 38 
mm and 51 mm tube diameters, the depth is assumed to have 
reached its ultimate value. From figure 3.12 to 3.16, it can 
be seen that as long as no slug forms, all of the average 
pressure data from different tube diameters fall onto the 
Same curve. The point where the graph deviates from this 
line indicates the beginning of the slug formation. This is 
verified from the flow visualization results. Hence the 
absence of the influence of tube diameters can be seen on 
the correlation equation developed. 

The present laboratory data was obtained using 
relatively small nozzle diameters. It might be expected that 
the characteristics of the small jets are different from 
those of a much larger jet (Davies et al., 1980). Therefore, 
any extrapolation to a larger «diameter jet for the 
correlation developed must be done in a cautious manner. 

Besides the basic nozzle design itself, the effect of 
the nozzle geometry was not investigated in the present work 
primarily because of a 90 degree elbow immediately upstream 
of the nozzle. This sharp turn creates secondary flow which 
reduces the coherence of the jet and is likely to nullify 


any gain in jet coherence due to careful nozzle design. 
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5.3 Differences between the simulated and actual down-hole 
geometry 

Many factors have not been incorporated into the 
present experimental model of the down-hole geometry. These 
include for example, orientation of the bedding planes, 
overburden pressure, porosity and permeability. Due to the 
inhomogeniety of the oil sands, complete modelling in the 
laboratory model would not in any event be possible. 

The formation of a slug of water in different tubes as 
in the case of 38 mm and 51 mm tube diameters depends on the 
volume flow rate and the standoff distance. Whether the slug 
will form in the real case can not be determined at this 
stage due to the unknown effects of the porosity and 
permeability of the oil sands. Only a very preliminary set 
of experiments on the effect of the porosity was performed. 
The permeability of this set-up was very low due to large 
hole size. Questions as to the size of holes and methods of 
obtaining representative values for porosity and 


permeability are by themselves of interest. 


5.4 Consequences of the predictions 

For a nozzle diameter of 6.4 mm and dynamic pressure of 
23 MPa, the ultimate depth of penetration is predicted to be 
5.9 meters and is achieved with a jetting time of 6 minutes. 
This result is very encouraging and shows the capability of 
using water jets to punch out holes in the oil sands 


formation. In addition to depths of penetration another 
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Quantity of interest is the water requirements. 

A plot of the water requirements versus dynamic 
pressure of the water jet to reach the ultimate depth of 
penetration can be seen in figure 5.5. The plot shows that 
as the dynamic pressure increases, more water is required to 
reach the ultimate depth of penetration. No optimum value 
for the water requirements is observed from this plot. It 
can be seen from this figure that there is no particular 
advantage in going to a much higher jet dynamic pressure for 
the water requirements will become excessive. For a dynamic 
pressure of 20 MPa, the ultimate depth of 5.4 meters is 
predicted with a jetting time of approximately 6 minutes. 
The water required to reach such depth is predicted to be 
2.3 m*?. For a large number of holes to be drilled by this 
water jet, the water requirements will certainly increase. 
Removal of this spent water will have to be considered in 
any future implementation of this type of project (eg. gas 
lift). A pulsed jet might be of great interest in reducing 
the water requirements. Chau (1980) found that a pulsed jet 
always give a higher penetrating power than a continuous jet 
and this is a very promising result. 

The model can not be conclusively evaluated because of 
the very limited field data available. No existing data on 
the ultimate depth of penetration on the actual tests can be 
found in any of the literature cited. To improve the 
evaluation of this model, a more complete set of field test 


data is needed. 
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Figure 5.5 An estimate of the water requirements to 
reach the ultimate depth of penetration for 
various dynamic pressures. 
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5.5 Alternative model in jet cutting 

The 'apparent' lack of fit for the low pressure results 
shown in figures 4.2 to 4.6 is due to the fact that the 
probability curves are being forced through one at 2/do 
equals zero. A better fit can be obtained by using a curve 


of the form: 


P(p>Pc)= C e Ml 2/de)* 


As can be seen, for example, in figure 5.6, at low 
pressures the probability of exceeding Pc is less than one. 
However, at some of the higher pressures this best fit gives 
values of the probability greater than one which is 
impossible. So then for the cases where C would exceed one, 
they are forced through one and hence the correlation is the 
Same as the original one for these cases. The values of A do 
not differ much between the first and second models as shown 
in figure 5.7. The influence of this second curve fit on the 
maximum depth of penetration is now depending on the values 
of the constant C and the probability assumed ® as shown 


below: 
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Solid wall (D=102 mm) 


do=3.2 mm A=/0.322x107* 
do=3.2 mm dA='0.274x10~-° 
Go=2.4 mm dA=.0.129x10-" 
Gdo=2.4 mm dA= 0.122x10-* 
do=1.6 mm A= 0.622x1075 
do=1.6 mm A= 0.480x107° 
Go=1t.6 mm d= 0.359x1075 
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Probability (Pressure > p critical) 
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The probability that the pressure at the piston 


ELUNE pic face exceeds the critical value for alternative model. 
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Figure 5.7 The variation of the parameter \ between the 


first and second model. ts 


5 ‘ a : : 
7 4 oa 

: : - i - _ o i 

ie. ! ; vt 

ie ‘ 4 vee" we 
7 area oP : 


j » i a 
Nir : 
a cs 1 
{ 
\ oy Sn) 
a A 
Tn, a 
: he a q ; 7 
) ue a: i nee oa i 
ARTI RMIT 
7 (a rie i ie 


i enh 
: 

6d SS 

ae | 


76 


For the case of a low pressure jet (2 MPa, do=3.2mm): 

Based on the first model: Zmax=765 mm 

Based on the second model : Zmax=739 mm (C=0.3) 
For the case of a high pressure jet (32 MPa): 

Based on the first model : Zmax=3514 mm 

Based on the second model : Zmax=3501 mm 

However the influence upon the cutting rate is more 
Significant particularly for the low pressure jets. If now 
the second curve fit to the probability is used, the cutting 


time becomes: 


It can be seen that the cutting times for the low 
pressure jets become substantially longer. For example, for 
a 2 MPa jet, C=0.3, so the predicted cutting time would be 
1.8 times longer than predicted by the first model. However 
for higher pressure cases the value of C is close to 1, e.g. 
for a 32 MPa jet, C=1.16 and the change in cutting time is 
14%. A summary of various different parameters for the first 


and second models can be seen in table 5.1. 
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Table 5.1 Variation of the parameters between the 


first and second model. 


Pa do D r Ms ‘ 


MPa mm mm xe 10) 


<s Cae 102 0.447 0.322 0.30 
oe Site TG2 0.394 0.274 hee: 
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0.5 16 102 0.046 0.048 1.06 
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6. Conclusions and Recommendations 


The major objective of this work was to determine an 
empirical correlation to estimate the maximum depth of 
penetration and the time to reach this depth for a high 
pressure water jet drilling into an oil sands formation. To 
this end a laboratory model was constructed to examine the 
characteristics of high pressure water jets jetting into a 
tube with one end closed. Flow visualization observations 
and records of the pressure at the bottom of the tube 
revealed that the water jet produced discrete impacts there. 
From this observation it was concluded that it is most 
appropriate to model the jet drilling as a stochastic 
process and to analyze the pressure signal statistically. By 
first requiring the probability of the pressure exceeding 
the critical value to be some non-zero value for 
non-negligible cutting rates to occur a correlation to 
predict the ultimate depth of penetration was developed. 
Secondly, by assuming a linear relation between ariiiing 
rate and the above probability, a second correlation 
predicting the depth of penetration with time was nevenoneds 
These correlations are: 


a. Maximum depth of penetration. 
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b. Depth at any instant in time. 


Dynamic pressure range: 15 MPa to 35 MPa. 
Nozzle diameter range: 1.6 mm to 3.2 mm. 


Standoff distance range: 0.3 m to 1.7 m. 


The results predicted from this correlation equation 
were found to agree fairly well with some limited field test 
data even though the assumed cutting model predicts the 
cutting rate to be independent of the pressure at the 
cutting face once the critical pressure is exceeded. A more 
complex cutting model was not pursued at this point for the 
laboratory data for establishing parameter values and field 
data for comparison purposes are both very limited. 

Extrapolation of the model to field scale conditions 
predict ultimate depths of 6 meters for a dynamic pressure 
of 23 MPa in 6 minutes which are certainly of interest. In 
addition, water requirements can be made with these 
correlations and a 2.4 m? volume of water is predicted to 
reach that ultimate depth. 

To further develop the stochastic approach utilized 
here to develop a cutting model, much more laboratory and 


field data of depth of penetration with time is required. 
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APPENDIX B 


The value of the parameter B was determined using the 
limited laboratory data from Chau (1980). Chau obtained 
depths of penetration under different jetting times (up to 4 
seconds) for various dynamic pressures. However, most of his 
data are limited to low dynamic pressures. From his 
observation, it was found that no cutting of oil sands 
occurs when the dynamic pressure falls below 1.2 MPa. Hence, 
to obtain a good estimate of the parameter B a higher 
dynamic pressure is essential. Therefore oe the depths of 
penetration for higher dynamic pressure were used in 
determining the parameter B. 


Using equation (9), that is: 


and the measured Gepth of penetration by Chau, different 
values of B were obtained for different jetting times and 
pressures. As a result, for each dynamic pressure several 
values of B were obtained, as shown in figure B.1. This 
means that whenever 1 second jetting time was used, the 
corresponding depth of penetration calculated using equation 
(9) for 1 second jetting time coincides with the measured 
depth of penetration. The same thing is true for 2 seconds 


etc. However, the calculated depth of penetration for 2 
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seconds using the value of B/obtained for 1 second will not 
give the same result as the measured depth observed in 2 
seconds shown in figure B.1. Therefore, the best values of B 
that will give the least error between the calculated depth 
of penetration and the measured values were determined. This 
was done by plotting the sum of the squares of the error, £, 
between measured depth and calculated depth to the values of 
B obtained for different jetting times as shown in figure 
Beene ine svalue of B tharewill give the minimum of this error 
summation was then chosen (ie. B4). 

It appears from the data shown in figure B.1, that this 
procedure results in a conservative (ie. low) estimate for 
the initial cutting rate. More extensive experimental data 


is required to improve this estimation. 
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Figure B.2 Determination of the -cuttingerate. 
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